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Abstract 

We study the fluctuations of a Brownian micro particle trapped with optical 
tweezers in a gelatin solution undergoing a fast local temperature quench below 
the sol-gel transition. Contrary to what was previously reported, we observe no 
anomalous fluctuations in the particle’s position that could be interpreted in terms 
of an effective temperature. A careful analysis with ensemble averages shows only 
equilibrium-like properties for the fluctuations, even though the system is clearly 
aging. We also provide a detailed discussion on possible artifacts that could have 
been interpreted as an effective temperature, such as the presence of a drift or a 
mixing in time and ensemble averages in data analysis. These considerations are of 
general interest when dealing with non-ergodic or non-stationary systems. 
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1 Introduction and Motivations 

1.1 Gelatin and the sol-gel transition 

Gelatin is a thermoreversible gel [1]. It is a heterogeneous mixture of water-soluble de¬ 
natured collagen protein chains, extracted by boiling animal by-products (skin, tendons, 
ligaments, bones, etc.) in water. Collagen molecules are rods of 300 nm length, made 
of three strands, with high average molecular weights. This triple-helix structure is sta¬ 
bilised by hydrogen bonds and has a diameter of ~ 1.4 nm. The chemical treatment used 
to produce gelatin breaks crosslinks between strands, but can also hydrolyze strands into 
fragments. Thus a broad molecular weight distribution is obtained for gelatin [2,3]. 

Above a temperature T^eit ~ 40 °C, the gelatin chains are in coil conformation. The 
gelatin solution is in a viscous liquid phase, called “sol” phase. Below a temperature 
Tgei ~ 30 °C renaturation of the native triple helix structure occurs, and chains form 
a percolating three-dimensional network of helical segments connected by single strand 
coils. The gelatin solution is in an arrested state with elastic behaviur, called “gel” phase. 
The coil-helix transition is completely reversible and the transition from one phase to the 
other is called the “sol-gel” transition [1]. 

Physical properties of the sol phase, and of the sol-gel transition are studied in [4,5] 
for different gelatin concentrations above 4wt %. In particular, it was seen that there 
are at least three successive steps in the transition: monomer to aggregate formation, 
random-coil-single-helix transition (disorder-order transition), and single-helix-triple-helix 
transition (order-order transition). It is then possible to identify different phase states in 
the sol domain: the sol state I where the chains have random coil conformations and the 
sol state II where single and triple helices begin to form (without reaching gelation). 

The gel phase was also shown to share properties with glassy materials, which are out- 
of-equilibrium metastable systems. After a quench at T < Tgei the system is frustrated by 
topological constraints because each gelatin chain is involved in at least two helices, and 
neighbouring helices are competing for the shared portions of non-helical chain. There¬ 
fore, the system displays physical ageing: its physical properties slowly evolve with time, 
through a process known as structural recovery. For example, the small-strain shear mod¬ 
ulus of a 5 wt % gelatin solution quenched at 20 °C increases logarithmically as a function 
of the ageing time [6]. And the elasticity of gelatin gels during slow cool and heat cycles 
exhibits memory and rejuvenation effects similar to the ones found in spin glasses [7]. 

Although it is known that mechanical properties of gelatin gels are very sensitive to 
temperature variations, previous thermal history of the gel, and time, this system has 


2 





some interesting experimental features: 

• The fact that the transition is thermoreversible allows us to do melting/gelation 
cycles simply by controlling the temperature of the sample. 

• The ageing rate can in theory be controlled by changing the quench depth. 

• The length-scale of the collagen chains (300 nm) is big enough to be sensed by a 
micro-particle of 2 pm. 

This particular sol-gel transition was chosen for previous works done in our group about 
fluctuations of Brownian particles in quenched gelatin samples [8-10]. 

A summary of the previous works results and our motivations are presented in the 
next section. 

1.2 Previous work: anomalous variance, heat flux and Fluctua¬ 
tion Dissipation Theorem violation in an ageing bath 

Previous works [8-10] showed that a particle trapped with optical tweezers in a liquid 
droplet of gelatin solution, quenched at a temperature below Tgei exhibits anomalously 
high position fluctuations right after the quench. These anomalous fluctuations could 
be interpreted as an effective temperature, and were consistent with a violation of the 
Fluctuation Dissipation Theorem (FDT) and with an exchange Fluctuation Theorem 
(xFT) for the heat exchanged between two heat baths at different temperatures [11]. 

We reproduce here some hgures from [8-10] and recall the associated key results: 

• The variance of the position = (x^) exhibits anomalously high value for short 
times (~ 5 s) right after the quench. It then stabilises at the equipartition value 
^Lq = ksT/k for 200 s (where ks is the Boltzmann constant and k the stiffness of 
the trap). And it hnally decreases logarithmically for long times after the quench. 
See hgure la. The anomalously high variances can be interpreted in terms of effective 
temperatures: a^. = /ceFeg/A: > cx^eq- 

• The Probability Distribution Functions of position fluctuations are Gaussian at any 
time after the quench, but their variances decreases with time (in agreement with 
the variances observed). See hgure lb. 

• The Probability Distribution Function of the heat exchanged between the particle 
and the bath during short times after the quench is asymmetrical. See hgure 2. 

• The asymmetry function of the heat p{q) = In (P(g)/P(—g)) (where P(g) is the 
probability of observing the value q of the heat) satishes an exchange Fluctuation 
Theorem: p{q) = A/3g, with its slope A/3 directly linked to the ehective tempera¬ 
tures dehned from the variances (see [9]). 

• The Fluctuation Dissipation Theorem is violated only for short times after the 
quench, and this violation can be linked with the amount of heat exchanged between 
the particle and the bath during the same time. See hgure 3. 

Unfortunately, none of those key results was found to be reproducible, and we believe 
that they were only due to an artifact in the data and/or in the analysis method. There¬ 
fore, we present in this article a detailed and careful analysis of trajectories of particles 
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Figure 1: (a) Evolution of the normalised variance ka'^/k^T of the position fluctuations 
of one particle trapped in gelatin solution (10 wt %) or glycreol, quenched at 26 °C, for 
different times t after the quench, (b) Evolution of the Probability Distribution Function 
of the position fluctuations of the particle trapped in gelatin solution. 





Figure 2: Probability Density Function of the normalised heat q exchanged during r = 30 s 
computed at different times t after the quench in gelatin. 


trapped in a droplet of gelatin solution quenched at a temperature below Tgei. We show 
that there is indeed no effective temperature for this system, which surprisingly exhibits 
equilibrium-like properties while aging. 

The article is organized as follow: in the hrst section we describe our experimental set¬ 
up, the gelatin solution preparation, optical trapping system, and local heating for fast 
quenching method. In the second section, we discuss our experimental results: we show 
with some bulk measurements that the system is indeed aging after the quench, we then 
analyze the effect of a slow drift in trajectory and of a mixing between time and ensemble 
averages in data analysis to show that there is no anomalous variance of the particle’s 
position. We end by discussing the consequences of the absence of effective temperature 
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Figure 3: Passive Power Spectral Densities of the position fluctuations (color points) 
and Fourier transform of the active response function (black dashed-lines) computed at 
different times after the quench, (a) For Os < t < 15s. (b) For 30s < t < 45s. (c) For 
75s < t < 90s. (d) For 1200s < t < 1215s. If the Fluctuations Dissipation Theorem 
is verihed, the two quantities should be equal, which is not the case for low-frequency in 
(a). 


for the heat exchange, Fluctuation Theorem and Fluctuation Dissipation violation. 


5 



























2 Experimental set-up 

2.1 Gelatin sample preparation 

We use gelatin powder from porcine skin, produced by Sigma-Aldrich®: gel strength 
~ 300 g Bloom, Type A, Bio Reagent, suitable for cell culture. This gelatin is derived 
from acid-cured tissue, whereas type B is derived from lime-cured tissue. 

We work with gelatin at a weight concentration of 5wt %. The samples are prepared 
following a standard protocol [12]: the wanted amount of powder is dissolved in bidistilled 
water, which is then heated for ~ 30 min at ~ 60 °C while slowly stirred until the solution 
is transparent and homogeneous. While the solution is still liquid, ~ 2mL are hltered 
using a Millex® syringe driven hlter unit with 0.45 pm pore size mixed cellulose esters 
membrane. Then 15 pL of an aqueous solution of silica beads (radius R = 1.00 ± 0.05 pm) 
with concentration 10^ particle mL“^ are added, and the solution is strongly agitated. The 
non-hltered and hnal solutions are let gel at room temperature and kept in the refrigerator 
for later use. 

We use a disk-shaped glass cell, with an Indium Tin Oxyde (ITO) coated microscope 
slide, a free-volume to avoid problems if the volume of solution changes during gelation, 
and a Wavelength Electronics TCS10K5 thermal sensor for temperature measurement (see 
hgure 4). To £11 the cell, the gelatin solution with dispersed silica particles is taken from 
the refrigerator and heated at ~ 50 °C until it is in the sol phase. 
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Figure 4: Schematic representation of the cell used to trap particles in gelatin solution 
(view from the side). The microscope slide is ITO-coated, which enables us to heat the 
cell by sending an electrical current through the glass surface. 


2.2 Optical trapping and controlled gelation 

For different purposes we used two different alignments of an optical tweezers set-up: 

The first set-up uses a laser beam (A = 532 nm) separated in two cross-polarized 
beams which enables us to have two traps with no interference between them. A laser 
diode (A = 980 nm) is aligned with the green laser and used to heat locally the sample. 
The tracking is done using a fast camera which is able to track two particles at 600 Hz. 
The acquisition rate is not very fast but by previously using a calibration target, we can 
directly convert the displacement of the particle from pixels to pm. See figure 5. 

The second set-up uses a single laser diode (A = 980 nm) to trap the particle and to 
heat locally the sample. A He-Ne laser (A = 632.8 nm) is aligned with the laser diode 
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Figure 5: Schematic representation of optical tweezers set-up used to trap two particles 
in gelatin solution with a green laser. An infrared laser diode is used to melt the gelatin. 
“M” is a mirror and “DM” are dichroic mirrors. 


and deflected by the particle. This deflection is measured using a position sensing diode 
which is able to track one particle at more than 10 kHz. Contrary to the first set-up, 
the position signal is in arbitrary unit and a supplementary calibration is needed for each 
measurement to convert the trajectory of the particle in physical units. See figure 6. 

For both set-ups, the microscope objective is an oil-immersion Leica HCX PL. APO 
x63 with high numerical aperture N.A. = 1.4. The microscope objective is surrounded 
with a custom-made heating ring, made with a Minco flexible resistor and a Wavelength 
Electronics TCS10K5 thermal sensor for temperature measurement. A feedback control 
is managed by the temperature control module (TCM-39032) of a modular laser diode 
controller (ILX Lightwave LDC-3900). As already mentioned, another thermal sensor 
is inserted directly inside the cell (see figure 4) and another feedback control is done by 
mastering the current going through the ITO-coated microscope slide, with an Instec MKl 
Board and a PID software. These two temperature devices ensure that the temperature 
of both the microscope objective and the cell are well controlled. The precision achieved 
on the temperature control is about ±0.05 °C. 

To trap particles, the temperature of both the microscope objective and the cell are 
set to 38 °C so that the gelatin is in the sol phase. Then, one particle is found and trapped 
at a given distance from the bottom surface of the cell (typically h = 15 pm). Then, the 
temperature controls are set to a value below the gel transition (typically Tfb = 27 °C) 
and we let the sample gel for a few hours (typically between 6 and 10 hours). Since the 
gelatin shows a lot of hysteretic behaviour [2,3], this gelation procedure appears to be 
important and the rheology of the gel can vary if the gelling time is very different {e.g. 
a few days). Moreover, one must pay attention to regularly check the distance between 
the bead and the bottom surface of the cell, since the focal distance of the microscope 
objective always drifts slowly when its temperature is changed. 

The refractive index of the liquid gelatin solution was measured Ugei = 1.3415, which 
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Figure 6: Schematic representation of optical tweezers set-up to trap one particle in gelatin 
solution. An infrared laser diode is used to trap and to melt the gelatin. The deflection 
of a He-Ne laser induced by the trapped bead is measured by the position sensing diode 
(PSD). The white light source and camera are only used for direct visualization but not 
for the measurements. “BS” is a beam-splitter, “IF” and interferential hlter to suppress 
the infrared beam and “DM” are dichroic mirrors. 


is close to value in water neater = 1.3335. It follows that the trapping stiffnesses in gelatin 
solution should be close to the ones in water with same experimental parameters. 

2.3 Local heating and fast quenching method 

When the particle is trapped and the sample is properly gelled at a given controlled 
temperature Tfb < Tgei, the local quenches are done in a similar way than the one presented 
in [8-10]: The power of the 980 nm laser diode is risen to a high value (typically^ 230 mW) 
during a given time (typically Tmeit = 200 s). Because of the light absorption of the water 
molecules in the solution, the temperature of the gelatin around the particle (which is at 
the focal point of the microscope objective) increases by a small amount 6T. Following 
the formula in reference [13]: 


6T 


Pa 

27iK 



( 1 ) 


^This is the power measured on the beam before the microscope objective, so the “real” power at the 
focal point should be smaller, due to the loss in the objective. 










































where a = 50 m ^ is the attenuation coefficient of water at 27 °C for wavelength 980 nm [14] , 
and K = 0.61Wni“^K“^ is the thermal conductivity of water. Here, we await: 

(5T~11°C. (2) 

This increase in temperature is only roughly estimated. Especially because we don’t really 
know what is the absorption of the microscope objective for the near infrared, and because 
it is impossible to measure the temperature with a usual probe on this very small scale. 
But it is seen that the increase is strong enough to melt a small droplet of gelatin (radius 
Ra ~ 10 pm) around the bead. Then, the power is quickly decreased to a low value (in 
the case where the same laser diode is used to trap and heat) or to zero (in the case 
where another laser is used to trap the particle), and the sample is let gel for a given time 
(typically Trest = 500 s). Since the thermal diffusivity of water is k = 0.143 x 10“®m^s“^ 
at 25 °C, the time needed to dissipate the heat from the droplet to the bulk is short: 

^ ^2xl0-^s. (3) 

Hence, the gelatin is believed to experience a fast quench at temperature Tfb < Tgei and 
should start ageing^. After the resting time Trest at low temperature Tfb, the power of the 
laser diode is risen again, and another quench is done. Note that the exact duration of 
Trest was uot Considered as important, because it was believed that the melting “resets” 
the gelatin sample and that all the anomalous behavior occurs right after the quench. 

The position of the particle trapped in the center of the melted droplet is continuously 
measured during a succession of several melting and aging. For each measurement the 
quenching is repeated a few hundred times in order to perform proper ensemble averages. 

An example of trajectory obtained with the second set-up is presented in hgure 7. 
When the intensity of the laser is high, the gelatin droplet is in the “sol” phase and the 
particle fluctuates in an optical trap with a high stiffness. When the intensity of the laser 
is low, the gelation is occurring and the stiffness of the trap is low (which is the reason 
why the position fluctuations are bigger). 


^Actually, in the case where the same laser diode is used for trapping and melting the droplet, the 
temperature of the quench is a little bit above Tfb because of the absorption of the laser. Since the 
power of the laser is low, this increase is less than 1 °C and can easily be compensated by lowering Tfb 
accordingly. 
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Figure 7: Trajectory of one particle trapped in a gelatin sample kept at Tfb = 27.5 °C, 
around the quench. On the red part of the trajectory, the intensity of the laser is high 
and the gelatin droplet is liquid. On the blue part, the intensity of the laser is low and 
the gelation is occuring. 


3 Results 

In this section we present some results showing that there are no anomalous fluctuations 
occurring in the ageing of gelatin solution right after a fast quench. We also discuss why 
this effect that was previously observed is likely to actually be an artefact due to data 
analysis. 

3.1 Time evolution of bulk properties and hysteresis 

We hrst performed preliminary measurement of gelatin gelation in bulk (he. without the 
local heating method). We prepare a cell with gelatin solution as described in 2.1. We 
set the temperature controls of both the objective and of the cell at 37 °C for 30 min to 
melt the gelatin. We trap one particle and we switch both temperature controls to a 
given temperature Tfb. We wait a few tens of minutes (typically ~ 30 min) and we set the 
distance between the particle and the bottom of the cell when there is no more drift due 
to thermal expansion of the microscope objective. After that, we measure the position 
of the trapped bead for a long time {e.g. 8h) at 400 Hz to see the bulk gelation of the 
sample. 

For all these measurements, the distance between the bead and the surface is h = 
15 pm. The liquid gelatin solution at 37 °C is a Newtonian fluid, and the stiffness k of the 
trap can be computed directly from the variance of the position x of the bead: 

= "f (4) 

For all these measurement the stiffness of the trap was k = 0.46 ± O.OlpN/pm. 
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Since the Tgei is expected to be around 29 °C, we varied Tfb from 31 °C to 27.5 °C. It 
was found that above 28.5 °C, the gelation does not occur on the time of the experiment 
and the solution stays liquid, even if its viscosity increases continuously. Below 28.3 °C, 
the gelation occurs before the end of the experiment. It was estimated that the bulk 
gelation of the cell volume takes ~ 260 min at 28.3 °C and ~ 120 min at 27.5 °C. 

Estimating the state (“sol” or “gel”) of the gelatin solution, is not trivial, since the 
fluid can be really viscous without being completely gelled. Qualitatively, the trajectory 
of the trapped bead starts to be heckled, and the bead sometimes escapes the trapping 
(see hgure 8 ). An a posteriori test consists in switching off the laser (resulting in switching 
off the trapping) and letting the sample at Tfb for a few more hours (typically over night) 
to see if the particle slowly fall to the bottom of the cell. If the particle does not fall, the 
gelatin solution is considered to be fully gelled. 


3 - 


2.5 


2 - 


o 

Oh 


1.5 


1 - 



200 210 220 230 240 250 260 270 
Time after quench (min) 


Figure 8 : Evolution of the position of one trapped particle, in gelatin solution kept at 
28.3 °C (after being melt at 37 °C for 20 min). At the end of the trajectory, gelation occurs 
and the particle is moved away from the optical trap. 


To estimate the evolution of the viscosity during the gelation process, we used passive 
micro-rheology techniques [15]. The trajectories were divided in portions of ~ Ih, and 
the Power Spectral Density (PSD) was computed for each portion. A long trajectory is 
required because we need low frequencies to correctly estimate the PSD. We explicitly 
assume that the bulk aging is slow enough for not perturbing too much the estimation 
of the PSD when a long trajectory is taken. Or at least, that taking a long time-window 
will only smooth the rheology result. 

As seen in hgure 9a, shortly after the decrease of temperature, the PSD is still 
Lorentzian, as awaited for a particle trapped in a Newtonian huid at equilibrium [16]. 
The viscous drag coefficient 7 = (with rj the dynamical viscosity of the solution) 

can be estimated from the value of the cut-off frequency fc = k/{2'K'y). Here we End: 
7 = 21 ± 1 X 10“^Pas. As the gelation occurs, the PSD is less and less Lorentzian (see 
hgure 9b), which is the sign that the gelatin solution starts to behave as a viscoelastic 
huid [17]. 

We plot in hgure 10 the evolution of the fitted cut-oh frequency at diherent time 
after the gelatin solution was set at Tfb = 28.5 °C. Even if the spectrum is no longer 
Lorentzian near the end of the measurement, it seems that the cut-oh frequency decreases 
exponentially. Therefore, the apparent viscosity increase is exponential. 


11 






(a) 70 min after the temperature change. (b) 320 min after the temperature change. 

Figure 9: Power Spectral Density of the position of one particle trapped in gelatin solution 
kept at 28.5 °C (after being melt at 37°C for 30 min). The PSD is estimated over a time- 
window of 1 h. Shortly after switching the temperature the gelatin solution is still a 
Newtonian fluid and the PSD is Lorentzian. After some time, visco-elastic effects appear 
and the PSD is no longer Lorentzian. 
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Figure 10: Evolution of the htted cut-off frequency fc, in gelatin solution kept at 28.5 °C 
(after being melt at 37°C for 30 min). 


From these preliminary measurements, we estimate that the Tgei is about 28.3 °C for 
our gelatin solution at 5wt %. We chose to work with Tfb < 28.3 °C for all the following 
quenching experiments. As mentioned earlier, gelatin solutions have big hysteretic be¬ 
haviour [2,3]. It follows that the viscoelastic properties of the solution in an important 
temperature range around Tgei cannot be known independently of the sample’s history. 

Another consequence of the hysteretic behaviour is that the hrst bulk gelation of the 
sample must be done in a controlled and reproducible manner. If the sample is let gel for 
a too long time (generally more than one day), or at a too low temperature (< 22°C), 
the hrst melting/regelling cycles used for the quenching experiment will be different from 
the following ones (where a “reproducible” state is reached). Especially, in this case the 
hrst melting is more difficult to reach. Examples of trajectories are shown in hgure 11. 
One can clearly see some position drifts occurring when the temperature is increased, 
before reaching a “sol” state where the particle huctuates in the optical trap. Note that. 
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(b) 3h after starting the cycles. 


Figure 11: Examples of trajectories for melting/regelling cycles with a gelatine sample 
gelled at Troom = 24 °C for a long time. The red dashed-lines indicate when the heating 
laser is switched ON, and the blue ones when it is switched OFF. At the beginning, the 
melting is more difficult to reach, after a given time, the cycles look “reproducible”. 


contrary to hgure 7, we use here the hrst set-up where the trapping laser is different from 
the heating laser (and the focal point of the two laser are not aligned). Hence, there is 
nearly no change of the trap stiffness when the intensity of the heating laser is changed. 
The change of fluctuations amplitude when the heating laser is switched OFF is mostly 
due to the rapidly increase of gelatin viscosity when gelation occurs. 


3.2 Difference between ensemble variance and temporal vari¬ 
ance in the presence of a drift 

We now consider quenching experiment as described in 2.3. We obtain several temporal 
trajectories of the particles positions for a given quenching temperature Tfb < Tgei. The 
important point is to estimate correctly the statistical properties from this set of data. In 
particular, we are interested in the variance of the position, which has been seen to have 
an anomalous increase right after the quench [8-10]. 

The correct ensemble variance should be estimated instantaneously at a given time 
t, by considering the N different trajectories at this time t. If one wants to increase the 
statistics by taking a small time-window 6t, there are at least 3 ways to compute the 
variance from the set of trajectories. These different ways are schematically represented 
on hgure 12. We call Xi{t) the position of the particle for the quench at the time t: 


• The temporal variance is obtained by estimating the variance over the time 5t 
for each quench, and then averaging over the N quenches: 


1 ^ 
i=l 


6t 


rt+St 

- Xi{t)) dt 


(5) 


where Xi{t) = ^ Xi{t') dt' is the temporal mean of x for the quench, between 
t and t + St. 


• The ensemble variance crg^gg^bie is obtained by estimating the variance over the N 
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quenches at a time t and then averaging over the time-window 6t: 


1 rt+St 


(J. 


ensemble 




N 




N-1 


i=l 


dt' 


( 6 ) 


where {x{t')) = is the ensemble mean of the N trajectories Xi{t') at 

time t'. 

The boxed variance is obtained by taking the N segments of trajectory from 
Xi(t) to Xi(t -f 6t), and then estimating the variance of the whole set of data: 


1 -'V 


(7) 


where [a[](t) = Xi{t') dt' is the mean computed on the set of data made 

of the N segments from Xi{t) to Xi{t + 5t). It is the variance used in references [8-10]. 

Nota Bene: Here to clearly distinguish the role of the time and the ensemble averages 
we have considered the time as a continuous variable and the number of trajectories as 
discrete. But experimentally the time is of course also a discrete variable, since we take 
measurements with a hnite sampling frequency. 
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Figure 12: Schematic representation of the different ways to estimate the variance for a 
set of N trajectories with a time-window St. The temporal variance is computed 
by estimating the variance of the points in the fuschia box, and then averaging over the 
trajectories. The ensemble variance is computed by estimating the variance of the 

points in the green box, and then averaging over the time-window St. The boxed variance 
(T^ox is computed directly by estimating the variance of all the points in the orange box. 


If the system is at equilibrium and St is big enough to correctly take account of the 
low-frequency of the motion, all these values should be equal to the equipartition value 
k^T/k, with /cb the Boltzmann constant, T the temperature and k the trap’s stiffness. 
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Unfortunately, when the system is non-stationary (which is the case for an ageing 
system), these 3 dehnitions of the variance are not equivalent. Especially, if there’s a slow 
drift existing on each trajectory, the estimations that average over time (he. temporal 
and boxed variances) are likely to show a strong artefact. 

To illustrate this effect, we have taken a set of 178 quenches done with the hrst set¬ 
up described in 2.2 at 28 °C, sampled at 400 Hz. The parameters were: melting time 
Tmeit = 250 s, melting intensity Jmeit = 235 W, resting time Trest = 305 s and trap stiffness^ 
k = 3.7pN/pm. One can clearly see on the trajectories that there is a small drift of 
~ 40 nm which occurs right after the quench (see hgure 13). Such a drift is often seen 
for this kind of measurement. We interpret it as a slow relaxation of the gel network, 
which occurs on a time much smaller than the gelation, but much greater than the heat 
dissipation. In other words, when the gelation occurs, the particle is trapped in the gel 
network at a given position. And even if we melt a small droplet, the gelatin network 
will somehow “remind” this position and pull the particle back to its place when it regels. 
Here the drift is very visible because the position of the trapping laser is not the same 
as the position of the locally heating laser. Thus the position where the particle was 
during the hrst bulk gelation is not the position where the particle is attracted to when 
the gelatin is melted. But even when there is only one laser used for both trapping and 
heating, this drift can occur. It is indeed impossible to verify that the position where the 
particle is when the sample gelled is exactly the position of the laser, and a drift of only 
a few nm can be visible. This kind of drift can be avoided by having a more powerful 
heating laser to completely melt the gelatin on a larger area, as in [18]. 



Figure 13: 20 hrst trajectories for a quench at Tfb = 28 °C, sampled at 400 Hz. A slow 
drift of ~ 40nm is clearly visible during the hrst ~ Is. After that, the position only 
oscillates randomly around a mean value. 


We then have three characteristic times : 

• Tgei the time needed for the gelatin solution to regel completely. It goes from a few 
hundreds to more than 1000 s depending on the quench temperature Tfb. 

• Aiynamics the typical time of the particle motion, which is directly j- and evolves 
from ~ 5 to ~ 100 s during the gelation process. 

^The trap stiffness is measured when the gelatin sample is completely melt and kept at constant 
temperature T = 37 °C, before the first bulk gelation. 


15 








• Tdrift the time where the drift is visible, which is typically Is for our experiment. 

If we take a 6t sufficiently small compared to Tdrift, the boxed and ensemble variances 
will give more or less the same result. Whereas, since Tdrift < '^dynamics, h is clear that 
the temporal variance will dramatically underestimate the variance due to the lack of 
low frequencies signal. Indeed, the temporal variance would require a 6t of the order of 
magnitude of Tdynamics for a correct estimation, which cannot be used because of the drift 
and the aging. Data are shown on hgure 14 for 6t = 0.1s. 



(a) Ensemble and boxed variances. 



Figure 14: Different variances computed for = 0.1s and normalised by the equilibrium 
value k-^T/k. The ensemble and boxed values are nearly equal and seem to be close to 
the equilibrium value at any time after the quench. Whereas the temporal value is clearly 
below the equilibrium value and decreases logarithmically with time after the quench. 


Now, if the chosen 6t is too big compared to the characteristic time of the drift, the 
boxed variances will start to show an anomalous increase. Data are shown on hgure 15 
for 6t = Is. This increase is not a real non-equilibrium effect due to the sol-gel transition, 
but only an artefact due to data analysis in presence of a slow drift. However, this slow 
drift is due to the fact that the sample is a gelatin solution, where an elastic network is 
created in the “gel” phase. 

It is nevertheless interesting to see that the correct ensemble variance seems to satisfy 
the equilibrium equipartition relation at any time after the fast quench, even though there 
is a clear evolution of the visco-elastic properties with time, and even in the presence of 
a slow drift at the beginning of the quench: 

Vt : t^ensemble(^) = (8) 

Similar results were seen for different quenches temperatures from 28 °C to 26 °C. 

In [8-10] it is stated that for longer times after the quench, the variance should decrease 
because of the elasticity of the gelatin network (as seen hgure la). This result is not clear. 
Of course, the particle dynamics becomes arrested after gelation, and at some point the 
amplitude of its positional huctuations must decrease in time. At the gel point, the storage 
modulus and the yield stress of the gelatin sample becomes non-zero, and even for a freely 
suspended particle the dynamics becomes subdiffusive [19]. Therefore it is clear that the 
temporal variance calculated for each trajectory should go to zero at long time. However, 
the ensemble variance calculated instantaneously over several trajectories should not go 
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Figure 15: Ensemble and boxed variances computed for a 6t = Is and normalised by 
the equilibrium value k-^T/k. The boxed variance clearly shows an anomalous increase at 
small times after the quench, which is the effect reported in previous works. 


to zero. Otherwise it would mean that the position where the particle stops to move is 
the same for every trajectory. Thus, the fact that the variance decreases for long time in 
previous works is a sign that the ensemble analysis is mixed with some temporal analysis. 

It was verihed experimentally that for Trest up to 900 s the correct ensemble variance 
remains constant. Afterwards, it is not clear that equation 8 is verihed, but we do not 
observe a shrinkage of the distribution of positions. 

3.3 Correct Estimation of Position Distribution Function 

We now want to study not only the variance of the position huctuations but also their 
complete Probability Distribution Function (PDF). In order to minimize the risk of slow 
drifts and to increase the sample frequency, we did new measurements with the second 
experimental set-up described in 2.2. With this set-up, the trapping and heating laser 
are the same, and the sample frequency can go up to 10 kHz. However, the calibration of 
the measured dehection from V to pm requires a supplementary assumption (for example, 
that the Fluctuation-Dissipation Theorem is verihed, as it is done in [20]). We will start 
by doing no assumption and plot the results only in arbitrary units. 

One must pay attention at which “position huctuation” is considered, as people often 
look at the distribution oi 5x = x — (x). The (x) is the mean of the position x, which 
can be dehned in several ways when the system is not a classic stationary ergodic system. 
Especially, when one considers a small time-window 6t, the correct mean should be the 
ensemble average {x{t)) estimated for each time t (as dehned in equation 6). But if one 
takes instead the temporal average Xi{t), estimated for each trajectory between t and 6t 
(as dehned in equation 5), the results will diher. 

As an example, we take the data of 132 quenches at 27.5 °C, sampled at 8 kHz. The 
parameters are: melting time Tmeit = 200 s, melting intensity Jmeit = 245 W, resting 
time Trest = 570 s, and trapping intensity /trap= 26 W which corresponds to trap stihness^ 
k ~ 5pN/pm. We compare the Probability Distribution Function of the positions with 
a 6t = 0.5 s at diherent times after the quench, when we subtract either the ensemble 

^The trap stiffness was measured in water (where viscosity is known) for the same laser intensity. 
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Position (a.u.) Position (a.u.) 

(a) When one subtracts the ensemble average, (b) When one subtracts the temporal average. 

Figure 16: Evolution of the Probability Distribution Function of the position fluctuation 
5x depending on the definition taken for the subtracted average. The PDFs are computed 
on a time-window 5t = 0.5s for different times after the quench, going from f = Os (blue 
curves) to t = 540 s (red curves). 


average (hgure 16a) or the temporal average (hgure 16b). In the hrst case, the PDFs 
are nearly always Gaussian and do not evolve in time. In the second case, the PDFs are 
always nice gaussians, but with a variance that decreases in time. This is consistent with 
the previous results showing that the ensemble variance is constant at any time after the 
quench, whereas the temporal variance decreases logarithmically with the time after the 
quench. And the variances estimated by doing a Gaussian £t on the PDFs clearly shows 
the same behaviour (see figure 17). 

This effect is simple to understand: the trajectories evolve on a time Tgei. This time is 
much bigger than tauc, the typical time of the fluctuations, and St. On the time-window 
St, each portion of trajectory Xi{t) can be written Xi{t) = Xi + Sxi{t), where Xi is the time 
average of the trajectory over the time-window. When one considers the N trajectory 
fragments between t and t + St, the difference between them is mostly due to the averaged 
value Xi of each trajectory fragment, and not to the fast fluctuations Sxi{t). Which means 
that the distribution of all the Xi{t) between t and t +St is nearly the same as the ensemble 
distribution of the x*. Whereas, the distribution of all the Sxi{t) is nothing more than 
the distribution of the fast temporal fluctuations of one single trajectory. 

This difference is very important, as any kind of high-pass hltering (for example a 
“detrend” function which is often used to suppress slow drifts) done to the trajectories 
will result in subtracting the temporal average, and thus distort the PDFs estimation. 

The experimental results show that the correct estimated PDFs do not evolve in time 
after a fast quench. Since we have already shown that the correct ensemble variance always 
verihes the equilibrium equipartition relation, we can conclude that the variance of these 
PDFs is simply k-gT/k. It is again interesting to see that, even if the gelatin solution is 
aging, its ensemble statistical properties seem to verify relations that are normally verified 
at equilibrium. 

It was also verified with some available data from [8-10] that the correct ensemble 
PDFs are not evolving with time after the quench. 
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Figure 17: Evolution of the variance estimated by fitting the PDFs with a Gaussian, at 
different times after the fast quench. When subtracting the correct ensemble average 
the variance is constant. When subtracting the temporal average, the variance decreases 
almost logarithmically with the time after the quench. 


4 What about heat and Fluctuation Dissipation The¬ 
orem? 

In previous works [8-10] the anomalous fluctuations observed right after the quench were 
interpreted in terms of heat exchanges between the bath and the particle. Indeed, the 
heat exchanged between t and f + r is equal to the variation of the particle’s energy 
/\Ut,r = AUt+r - AUt: 


Qt,r = AUt,r = ^(x‘^{t + T)-x‘^{t)y (9) 

In particular, the decrease of the variance after the quench has been interpreted as the 
sign of a heat transfer from the particle to the bath : 

(<5t,r) = ^ (c^^(t + r)-a^(f)) < 0. (10) 

The Probability Distribution Functions (PDF) of the Qt^r were shown to be asymmetrical 
for values of t and r chosen right after the quench (he. where the anomalous fluctuations 
were observed). 

A violation of Fluctuation-Dissipation Relation was also observed for times right after 
the fast quench. It was linked to the non-zero heat exchange by a modihcation of the 
Harada-Sasa equality [21,22] for non-stationary systems: 



^^Im{R(f,/)} 


d/ 


2|(gt,At)i 

k 


(11) 


Where Sx(t, f) is the Power Spectral Density of x and R(t, f) is the Fourier transform 
of the linear response function of the position a: to a perturbative time-dependent force 
(these two quantifies are function of the frequency /, but also of the time t since the 
system is ageing). 
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All these interpretations conies from the fact that the variance was seen anomalously 
high right after the quench, and then reduces to the equipartition value after a given time. 
In particular, the asymmetry and the shape of the PDFs of Qt,T are simply mathematical 
consequences of the fact that x{t+T) and x{t) have Gaussian PDFs with different variances 
+ r) > cr^(t). For example, the exchange Fluctuation Theorem (xFT) that was 
retrieved with the asymmetrical PDFs of Q is mathematically verihed for any random 
variable dehned hy y = xi — X 2 where xi and X 2 are random variables with centered 
Gaussian distribution of different variances: ^ 

Since we have already shown that, if estimated correctly, the PDFs of x show no 
anomalous behavior and have a constant variance equal to k-^T/k, if follows directly that 
the PDFs of Qt,T are symmetrical. Gonsequently, in average there is no heat exchange 
between the particle and the bath, for any t and t + t, and the xFT reduces to a trivial 
equality because the asymmetry function of a symmetrical distribution is always zero. 

Gonsidering the Fluctuation-Dissipation Theorem, one must remind that it is a priori 
not a good idea to test it in Fourier space. Indeed it is necessary to assume that the 
system is stationary and ergodic to link the correlation function to the power spectrum 
with the Wiener-Khinchin theorem [23,24]. Therefore, when the system is not stationary, 
one should in theory look at the proper ensemble correlation function: 

1 ^ 

EnsGorr^ 3 .(t,r) = [xiit) - {x{t))] x [xi{t + r) - {x(t + r))] (12) 

instead of the Power Spectral Density (PSD), which is a temporal quantity. Of course, 
one can always dehne a PSD of Xi on a given time-window 5t for each trajectory S'j;.(f, /). 
And this PSD would be equal to the Fourier Transform (FT) of the temporal correlation 
of Xi computed on the same time-window: 

1 rt+St 

TimeGorra;a;(f, r) = — [xi{t') - Xi] x [xi{t' + r) - Xi] dt' = FT{Sxi{t, /)}. (13) 

But the system needs to be considered stationary and ergodic on the time-window 6t, so 
that the ensemble and temporal correlations should be equal. 

Here, the assumption of local stationarity is reasonable since the PSD were computed 
on 15s long time-windows (which is short compared to the ~ 900s necessary to gel). 
However, it seems probable that the observed violation of Fluctuation-Dissipation Theo¬ 
rem was only due to the same kind of artifact already responsible for anomalous variance 
increase (for example: slow drifts for times right after the fast quench), because PSDs are 
sensible to low-frequency noises. Thus, there is no reason that this apparent violation is 
linked to an heat exchange, which does not exist anyway. 

Some ensemble correlation functions of the particle’s position are shown in figure 18 
for a set of 40 quenches at 27 °G, sampled at 8 kHz. The parameters are: melting time 
T'meit = 200 s, melting intensity = 270 W, resting time Xrest = 570 s, and trapping 
intensity /trap= 26 W which corresponds to trap stiffness^ k ~ 5pN/pm. The data are very 
noisy, but there is a tendency: the characteristic time increases after the quench, which is 
reasonable since the gelatin viscosity is also increasing during the gelation. The correlation 
functions are not simply exponential relaxations, which is consistent with the fact that 
the PSD are not Lorentzian (as shown in hgure 9b). 

We also made some experimental tests of Fluctuation Dissipation Theorem (FDT), by 
looking at the ensemble correlation of x and the response to an Heaviside change of the 

®The trap stiffness was measured in water (where viscosity is known) for the same laser intensity. 
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Figure 18: Normalised ensemble correlation function for a quench at 27 °C. Here we 
keep t fixed and we vary r from —10s to Os. The normalisation is done by dividing 
EnsCorr 3 ; 3 ;(t, r) by the value of k-^T/k extracted from the variance of the position PDFs. 


position of the trap. For these measurements, the position of the trap is changed from Xi 
to X 2 at a time Ir after the hrst quench, and the sample is let gel with the particle in X 2 . 
Then, for the second quench, the position of the trap is moved back to Xi at time Ir after 
the quench, and the sample is let gel in Xi. The procedure is then repeated alternatively. 
The perturbation introduced by the change of trapping position allows us to compute a 
normalized response function, averaged over the trajectories: 


X{tRX) = 


{x{tR + r) - ^initial) 


X, 


final 


X 


initial 


(14) 


where [Xnit.iai; X finai ] = [Xi;X 2 ] or [X 2 ;Xi]. It corresponds to the usual dehnition of the 
response function: 

_ ( 3 ^(^)perturbed 2:(t)unperturbed) (15) 

perturbation amplitude 

We use Xfinai—Xinitiai which is proportional to the perturbation amplitude. And we simply 
take Xinitiai as the average value of the unperturbed trajectory, because the mean position 
of the bead is constant and equal to the position of the trap if there is no perturbation®. 
If the FDT is verihed, the response function should verify: 

xitR, r) = 1 - jX^EnsCoTT^^itR, r) (16) 

kuT 


Some data are presented in hgure 19 for 50 quenches at 26 °C, sampled at 8 kHz. The 
parameters are: melting time Tmeit = 200 s, melting intensity Jmeit = 270 W, resting 
time Trest = 570 s, and trapping intensity /trap= 26 W which corresponds to trap stiffness^ 
k ~ 5pN/pm. The values of Xi and X 2 are estimated by computing the mean position 
of the bead when the gelatin is melted (which gives alternatively Xi and X 2 ). The exact 
value of k-^T/k was extracted from the variance of the position PDFs computed before 
changing the position of the trap. These measurements are a bit noisy because it requires 
a lot of statistics to compute a proper ensemble correlation function, but no apparent 
violation of the FDT was found for the times tested. 

®One could also take {x{tR)) to guarantee that x(tK,0) = 0, but it wasn’t necessary here. 

^The trap stiffness was measured in water (where viscosity is known) for the same laser intensity. 
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100 s after the quench 



Figure 19: Normalised response function xi^Ri '^) ensemble correlation function for 
tR = 100 s after the quench, and r going from 0 to 10 s. 


We didn’t test the Fluctuation Dissipation Theorem for times taken shortly after 
the quench, because the ensemble correlation shows a characteristic time which is very 
short at this time (see hgure 18). It is then more difficult to compute a proper ensemble 
correlation right after the quench, than when the viscosity of gelatin has already started 
to increase. We also didn’t compute the response function by varying tR for a hxed t_R + r, 
because it would require a lot of time to do the experiments. Indeed, each set of tR requires 
one day of measurement to compute r), and the sample cannot be kept a lot of days 
without degrading. 
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5 Conclusion 


In conclusion, we have locally studied the gel transition of a gelatin solution. We were 
unable to reproduce the results of previous works [8-10], but we have identihed some ex¬ 
perimental and data analysis artifacts which may explain the effects previously observed. 
In particular we have analyzed the effect of time-windows on proper ensemble averages, 
which are important to study aging systems. 

We have shown that in the hysteresis range of temperature (28.3 °C < T < 36 °C), 
bulk gelation can occur on very long times, and viscoelastic properties gradually appear. 
The characteristic time of the particle trapped in the bulk-gelling sample was seen to 
decrease exponentially before the gelation (whereas the viscosity evolves logarithmically 
after the gelation). 

For fast quenches of a small droplet of gelatin solution, we have found that the Prob¬ 
ability Distribution Functions of the position of the trapped particle do not evolve with 
time after the quench, even if the gelatin sample is undergoing aging and the viscoelastic 
properties are clearly evolving. Moreover, these PDFs show equilibrium-like properties, 
being Gaussian with a variance equal to the equipartition value k^T/k. These results 
seem not so surprising a posteriori, since it was already observed in the previous works 
that, after ~ 15 s the Brownian motion of the trapped particle behaves like in equilibrium 
with the thermal motion of the gelatin chains. Only the very first seconds after the quench 
showed anomalous behavior, which was strange, because the complete gelation occurs on 
much larger scales (~ 900 s). It however remains striking that a system which is clearly 
not stationary because of aging has ensemble properties which are stationary. 

For systems which are not ergodic or not stationary, time properties can be very 
different from ensemble properties. And it was also shown that some artifacts (like slow 
drifts) or analysis bias (like high-pass filter) can greatly modify the results if ensemble 
properties are estimated on time-windows. Therefore, one must be very careful when 
studying statistical properties of an aging system. This kind of problems had already 
arisen for other aging systems. For example, it was already shown in [25] that increase in 
effective temperature previously seen in suspension of Laponite [26] were in fact artifacts 
due to analysis methods. 

Finally, in agreement with the absence of anomalous position fluctuation after the fast 
quench, no heat exchange, nor violation of the Fluctuation Dissipation Theorem was seen, 
as it would be expected in an equilibrium medium. 
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